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ABSTRACT 

The possible existence of (meta-) stable de Sitter vacua in string theory is of 
fundamental importance. So far, there are no fully stable solutions where all 
effects are under perturbative control. In this paper we investigate the presence 
of stable de Sitter vacua in type II string theory with non-geometric fluxes. 
We introduce a systematic method for solving the equations of motion at the 
origin of moduli space, by expressing the fluxes in terms of the supersymmetry 
breaking parameters. As a particular example, we revisit the geometric type 
II A compactifications, and argue that non-geometric fluxes are necessary to 
have (isotropically) stable de Sitter solutions. We also analyse a class of type 
II compactifications with non-geometric fluxes, and study the distribution of 
(isotropically) stable de Sitter points in the parameter space. We do this 
through a random scan as well as through a complementary analysis of two- 
dimensional slices of the parameter space. We find that the (isotropically) 
stable de Sitter vacua are surprisingly rare, and organise themselves into thin 
sheets at small values of the cosmological constant. 
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1 Introduction 

Since the turn of the millennium, there have been many attempts to embed dark energy and 
inflation in string theory using flux compactiflcations. For this purpose, a whole class of string 
theory constructions involving flux backgrounds compatible with minimal supersymmctry 
has been considered in the literature. In particular, the mechanism of inducing effective 
superpotentials from fluxes [1] has often been considered in the context of the supergravity 
descriptions normally referred to as S'TfZ-models [2-10]. It has, however, turned out to 
be surprisingly difficult to find the required stable or quasi-stable de Sitter (dS) solutions 
through compactiflcations of string theory. Let us start with a brief review of the main 
results of these searches. 

Based on the results of [11], in ref. [12] all the type llA coset compactiflcations with SU(3) 
structure were analysed, and a no-go result was proven in all the cases but one. The only 
model evading these no-go theorems turned out to be an SU(2) x SU(2) compactiflcation. 
This model was indeed numerically shown to contain a dS solution, which was found to 
be unstable due to the presence of a tachyon already within the isotropic subsector. This 
example of a type llA unstable dS extremum was later studied in ref. [13] (see also [14]), 
where it was shown to belong to a whole continuous line of tachyonic solutions, for which the 
ten-dimensional uplift to massive IIA supergravity was constructed in terms of very simple 
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torsion classes. This was accommodated in a general framework in [15]. For further analysis 
of this model, see ref. [16]. In parallel, in ref. [17], many type II setups in the Z2 x Z2 orbifold 
with generalised fluxes were analysed searching for dS solutions. In a setup corresponding to 
type IIA with 06/D6 and only geometric fluxes, a dS region was found in the neighbourhood 
of a Minkowski critical point, but again with no room for stability. Later in ref. [18], the 
original dS solution of ref. [12] was localised at a particular point within the aforementioned 
dS region. 

Finally in ref. [19], the complete set of solutions in the framework of isotropic J\f — 1 
geometric type IIA compactiflcations on T^/ (Z2 x Z2) was presented. It turns out to consist 
of a discrete set of hues, of which only one crosses the Minkowski point [20] and enters the 
dS region, with tachyons at every point. This can be regarded as an indication that one 
needs to consider more general setups including non-geometric fluxes [21] in order to flnd 
stable dS vacua [17] . 

The main goal of this paper is to show how non-geometric fluxes indeed provide enough 
freedom in parameter space for tuning the masses of all the moduli to be positive. Moreover, 
in situations where this is the case, we would like to locate the region of stable dS in parameter 
space and qualitatively estimate which portion of the total volume this corresponds to. 

The paper is organised as follows. In section 2 we introduce the Z2 x Z2 compactiflcations 
with orientifold and their superpotential formulations. Subsequently, we present a new 
technique for solving the equations of motion in the origin of moduli space. In such a 
framework it becomes straightforward to predict the dimensionality of the parameter space 
of solutions. Finally, we make use of this machinery to review type IIA compactiflcations 
with geometric fluxes regarded as an over-constrained particular sub-case. In section 3 a 
IIB setup with non-geometric fluxes is introduced in which the parameter space has the 
correct degrees of freedom to gain control of the masses for all the isotropic moduh. A 
statistical scan of the parameter space is performed, and the fraction of stable solutions is 
plotted as a function of the normalised cosmological constant. Furthermore, we show some 
two-dimensional slices of the parameter space, and locate the thin region of stable dS critical 
points. Finally, we present our conclusions. An appendix is provided where we give some 
details about the analytical solutions in type IIB with F, H and Q fluxes. 
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2 Af = 1 compactifications with generalised fluxes 



2.1 The Z2 X Z2 orbifold 

A wide class of string compactifications down to four dimensions compatible with minimal 

supersymmctry can be effectively described by particular supergravity theories which are 
known as STU-mode\s. These theories arise, e.g., from j{%2 x ^2) orbifold compactifica- 
tions of type IIB with 03/07-planes (and duals thereof). This class of string compactifica- 
tions has the very appealing feature of collecting together low-energy effective descriptions 
which are related to each other by dualities without the intervention of mirror symmetry. 
In the Z2 X Z2 orbifold any dual effective theory is still formally described by the same 
Lagrangian where only fields and couplings have been transformed. 

This class of effective descriptions enjoys M = \ supersymmctry and SL(2)'' global 
bosonic symmetry due to the coupling of the gravity multiplet with seven extra chiral multi- 
plets. The scalar sector contains seven complex fields spanning the coset space (SL(2)/SO(2))''' 
which we denote by = (S*, Tj, with i = 1, 2, 3. In the type IIB language, the S mod- 
ulus contains the ten-dimensional dilaton, whereas the T and U moduli are interpreted as 
Kahler and complex structure moduli respectively^. The kinetic Lagrangian can be derived 
from the following Kahler potential 

3 3 
X = - log {-i {S -S)) -Y, log (-i (T, - T,)) - log (-i {IJ, - Ui)) . (2.1) 

i=l i=l 

This yields 

dSdS A / OTidTi dUidUi \ , , 

Ain = —2 + Y 7 - ,2 + 7 - ,2 • 2.2 

The presence of fiuxes induces a scalar potential V for the moduli fields which is given 
in terms of the above Kahler potential and a holomorphic superpotential W by 



V 



= (^-3\W\^ + K""'^ DaW D-pW^ , (2.3) 



where is the inverse Kahler metric and D denotes the Kahler-covariant derivative. 

As a first example of fiux-induced deformation, let us consider the presence of gauge fiuxes 
in the duality frame of type IIB with 03 and 07-planes which gives rise to the following 
superpotential 

W = Pf{U,) + PnmS , (2.4) 

F flux H flux 



^Please note that, in a type IIA language, the ten-dimensional dilaton would correspond to a combination 
of S and T, whereas the Kahler and complex structure moduli would be interchanged. 



where Pp and Ph are cubic polynomials in the complex structure moduli given by 

i i 
i i 

The Af — 1 supergravity defined by the above superpotential has a no-scale feature due to 
the absence of the moduli Tj. This implies that their real parts appear as completely fiat 
directions in the scalar potential. The line of including some non-perturbative effects such as 
gaugino condensation [22] has been considered as a possible mechanism to further stabilise 
the Kahler moduh (see e.g. refs [23,24]). 

Duality covariance arguments would suggest to complete the superpotential given in (2.4) 
with more effective couplings which are obtained by performing SL(2)^ transformations. This 
gives rise to a generalised form of W containing all the possible 2^ = 128 couphngs from the 
constant one up to the ST1T2T3U1U2US term. Such a very general superpotential, though, 
will be non-geometric in any duality frame. 



2.5 



The isotropic truncation 

For the sake of simplicity, in the main part of this paper we will restrict ourselves to the 
isotropic limit of the theories previously presented. This reduces the number of complex 
moduli down to three via 

Ti = T2 = Ta = T and = U2 = Us = U , (2.6) 

which span the scalar coset ( gQ^2) ) ' fiuxes, taking the isotropic limit implies 

the following identifications 

a['^ = af^ = af^ = (2.7) 

and similarly for ag'', hf' and 62^. 

In the isotropic limit, the complete duality covariant superpotential contains 32 couplings 
going from the constant term up to the ST'^'U^ term. For completeness, we give here the full 
form of the flux-induced W and its stringy interpretation in the different duality frames. 

The orbifold involution forces the six-dimensional internal space to be factorised into three 
2-tori as shown in figure 1. Now we will summarise the correspondence between generalised 
isotropic fiuxes and superpotential couplings appearing in the J\f — 1 effective description. 
The generalised fiux-induced superpotential can be written as 

W={Pf + PhS) + 3T{Pq + PpS) + 3T^Pq, + Pp,S)+T^Pf, + Ph'S) , (2.8) 
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Figure 1: = x T| x T| torus factorisation and the coordinate basis. We will use indices 
a,b,c for horizontal "— " directions 1, 3, 5 andindicesi,j,k for vertical "|" directions 2, 4, 6 
in the 2-tori with i — 1,2,3. 

where^ 

Pp = ao-3aiU + 3a2U^ -asU^ , Ph = ho - 2>hiU + UiU"^ - hU'^ , 

(2.9) 

Pq = co + ciC/-C2C/2-C3C/3 , Pp^do + diU -d2U^ -dsU^ 

and 

Pp, = a^5 + 3a'2t/ + 3a;t/2 + a'o[/3 , Ph' = b'^ + Sb'^U + Sb'^U^ + b'^U^ , 

(2.10) 

Pq, = -c^ + 4 [/ + ci C/2 -c'oU^ , Pp' = -d^ + 4 [/ + < -d'^U^ . 

The first half of the couphngs represent fluxes which admit a locally geometric interpreta- 
tion in type IIB (unprimed fluxes), whereas the remaining ones (primed fluxes) represent 
additional generalised fluxes which do not even admit a locally geometric description. These 
were first formally introduced in ref. [9] as dual counterparts of the unprimed fluxes. 

For our present purposes we focus on superpotentials induced by generalised fluxes in 
type IIB with 03 and 07-plancs which describe a locally geometric string background. Such 
a form of the superpotential turns out to be given by 

^(loc.geom.) ^ (p^ ^ ^ 2,T {Pq + Pp S) , (2.11) 

where each polynomial Pp, Ph, Pq and Pp has degree three with respect to U . Their explicit 
form and their interpretation in terms of the IIB fluxes are respectively given in (2.9) and 
table 1. 

In this paper we will be interested in constraining the stability of solutions of these 
isotropic S'TfZ-models originating from (generalised) string compactifications. The first rel- 
evant constraint comes from the analysis of the supersymmetry breaking sector and its 

^Please note that, in principle, the truncation to the isotropic sector gives rise to 32 + 8 = 40 fluxes, where 
all the fluxes transforming in the mixed symmetry representations of GL(6) (i.e. Q, P and their primed 
counterparts in type IIB) have in fact two fluxes (ci,ci) etc. giving rise to one single coupling (2ci — Ci) 
etc., so without loss of generality, we set throughout the text ci — ci etc.. 
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Table 1: Mapping between unprimed fluxes and couplings in the superpotential both in type 
IIB with 03 and 01 and in type IIA with 06. The six internal directions depicted in figure 1 
are split into " — " labelled by i = 1,3,5 and " | " labelled by a = 2, 4, 6. Note that the empty 
boxes in type IIA are related to the presence of dual fluxes analogous to the 'primed' notation 
in type IIB. The last column refers to the number of independent superpotential couplings 
the corresponding fluxes give rise to once the non-isotropic moduli are opened up. 

universal features which allow one to constrain metastable de Sitter vacua [25] and infla- 
tion [26,27]. The averaged projection of the scalar mass matrix along the supersymmetry 
breaking directions results in the following sGoldstino bound for the 7] parameter 

77 = Min jEigenvalues (^(m^)^^ j I < ^ " ' (2-12) 

where 7 = g^^^^p and IZ denotes the sectional curvature, i.e. the projection of the Riemann 
tensor along the sGoldstino. 
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Tabic 2: Mapping between primed fluxes and couplings in the superpotential. The conventions 
are the as in table 1 and again, just as there, the empty column should be filled in with extra 
dual fluxes. 

In the present supergravity setup we can explicitcly evaluate the range of variation of the 
sectional curvature term in (2.12). TZ can be written as 2/neff, where 

{GsGSf + \ {GrG^f + | {GuG^f ' ^ ' ' 

where G denotes the Kahler-invariant combination [28] K + log and its first derivative 
Ga is normalised such that GqC" = 3(1+7). ries gives a notion of how many fields effectively 
participate in the supersymmetry breaking phenomenon. Such a quantity turns out to be a 
real number between 1 and 7, these two extremal cases corresponding, respectively, to the 
limits in which either Ga only points in the 5" direction or Gq, ~ (1, 3, 3)^. 

■^The T and U components are equal to 3 instead of 1 because this model originates from the SL(2) theory 
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The sGoldstino bound in (2.12) implies in this case that, whenever V > 0, 



neff> 3(1 + 7), (2.14) 

which roughly tells us that stable dS favours homogeneous supcrsymmctry breaking rather 
than the other extremal situation in which it takes places mostly in one complex direction. 

2.2 A new systematic prescription 

Now we will present a general method for symplifying and subsequently solving the field 
equations in the origin of moduli space. If the scalar manifold is homogenous and the 
superpotential couplings constitute a closed set under non-compact duality transformations^, 
restricting the search for solutions to the origin of moduli space does not imply any loss of 
generahty. This argument was first employed in ref. [19] in the context oi J\f — A and J\f — 1 
supergravity. The advantage of this is that the equations of motion for the scalar fields turn 
into a system of quadratic conditions for the fluxes. 

Our prescription relies on a linear change of variables which introduces a number N of 
supersymmetry breaking parameters, N being the number of real fields in the model. The 
"broken" supersymmetry equations in the origin $o for the complex moduli {^'^}a=i n/2 
read 

DaW\^, = + iB^ , (2.15) 

where and denote N real constants and the Ihs of (2.15) is a linear expression in the 
M fluxes. This imphes that, whenever M > N, the above equations can always be solved in 
terms of fluxes as a linear function of the supersymmetry breaking parameters. 

By plugging the expression of the solution of (2.15) into the field equations in the origin, 
one obtains a new system of equations in which the new unknowns are the N supersymmetry 
breaking parameters and the remaining (M — N) fluxes. The peculiar feature of this new 
system lies in the fact that the extra fluxes will now appear only linearly in the fleld equations, 
the only unknows which still appear quadratically being A^, and Ba. If one arranges the 
supersymmetry breaking parameters into a real vector Xj with I — 1, . . . , N and denotes 
the fluxes by {F^} ^, the fleld equations will take the form 

h'^XiXj + q'aXiF^ = , (2.16) 

in the isotropic limit, where the triplets of complex moduli Tj and Ui collapse into T and U respectively. 

'^Please note that the homogeneity condition is satisfied in supergravity theories arising from Z2 x Z2 
orbifold compactifications, both in the isotropic and in the non-isotropic case. As for the second assumption, 
it turns out to be true in the two examples which are relevant here, i.e. type IIA with geometric fluxes and 
type IIB with generalised fluxes. 
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where h and g denote some real constants. The reason why this always has to be the case 
is that the supersymmetry equations should imply the field equations. Hence, whenever the 
rhs of (2.15) vanishes there cannot be any quadratic conditions left to be imposed on the 
remaining fiuxes. 

This makes it particularly simple to solve the field equations in the origin in the case in 
which the extra fluxes remaining after the introduction of the supersymmetry breaking pa- 
rameters are at least N (this requires M > 2N) . This sequence can be sketchily represented 
as 



M 



(2A5) 



\ 



?M 



\ 



?2N 



( F' \ 



p2N 



F^) 



J 



A situation which then turns out to be of special interest is that of having exactly 2N 
fluxes turned on, so that one is able to write down the general solution to the fleld equations 
by expressing all the fluxes as a function of the supersymmetry breaking parameters. The 
A^-paramater family of solutions generated in this way should provide, at least as a matter 
of principle, enough freedom to have control upon the eigenvalues of mass matrix for the N 
scalars retained in the theory. Generically, this tachyon-free region in parameter space will 
have a non-trivial intersection with the dS region. This situation is conceptually depicted in 
flgure 2. 

Our perspective points towards the possibility that, by introducing a non-isotropic setup 
admitting a 14-dimensional space of solutions, one might have enough freedom for giving 
a positive mass to all the scalars in the SL(2)^ model. At least qualitatively, increasing 
the number A^ of fields should not spoil this picture. The issue that we will address in the 
next section will be extimating how small the region of stable dS is compared to the whole 
parameter space of solutions. This will tell us how difficult it practically becomes to scan 
the whole parameter space searching for stable dS vacua. This analysis, on the contrary, 
will certainly depend on A^; in order to make predictions, we will now focus on the simpler 
N — 6 case corresponding with the isotropic model. 



2.3 Geometric IIA backgrounds: an over-constrained example 

Let us now consider the type IIA duality frame with 06-planes. This class of compactifica- 
tions fits again the framework of A/" = 1 effective descriptions in that there are two different 
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Figure 2: The dS region as well as the tachyon-free region are subsets of positive measure in 
the N -dimensional parameter space of solutions. Generically, they will overlap in a region 
that may be small but still is expected to be of positive measure. 

kinds of allowed D6-branes, i.e. those ones parallel to the 0-planes which do not further 
break supersymmetry and those ones orthogonal to them which break supersymmetry down 
to M' = 1. In this context, in the isotropic limit, one can have the four different R-R gauge 
fluxes ({-Fpjp^Q 2 4 6)' ^^'^ independent components of H flux, and three components of met- 
ric flux u, of which only two generate independent superpotential couplings (see footnote 2). 
The corresponding superpotential reads 

^(geom.iiA) _ {ao-3aiU + 3a2U^ -asU^) + {bo-3biU)S + 3T{co + ciU) , (2.17) 

where the type IIA interpretation of the above couplings can be found in table 1. In this case, 
setting ci = Ci as discussed in footnote 2 has a special geometric interpretation, because it 
represents the correct way of solving the constraints u"^ = and u H = which ensure the 
closure of the exterior derivative defined by the corresponding NS-NS background. 
We choose the origin of moduli space to be 

So = To = Uo = t . (2.18) 

By applying the counting introduced in section 2.2, one immediately realises that the amount 
M = 8 of fluxes is not large enough for following the above prescription for solving the equa- 
tions of motion in the origin. After introducing the six supersymmetry breaking parameters 
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by solving the (2.15), one finds 



02 

bo 
bi 

Co 
Cl 



1 (-Bi -B2 + 2Bs) 

I (-6^1 + 3^2 + A3) 
I (25i - ^2 - 53) 
-| (3^1 + A3) 



(2.19) 



= -|(Si + 53) 



After solving two axionic equations of motion by Oq = Oi = 0, one is left with the following 
field equations 



-3Al + 2A1A3 - 2A2A3 + Al- Wl - 2B1B3 - 2B2B3 + Bi 
9Al - 6Ai{A2 + As) - 2A2A3 + SAl + 9Bf - - 6S1S3 - 6S2S3 + 3S| 

9AI - 9A1A2 + SA1A3 + 3A2A3 + 3Bf - 3B1B2 - B,B3 - B2B3 - ABl 

6A1B1 - 3A1B2 + 3Ai53 - 3^251 + 3/12^3 + ^3^1 + A3S2 










(2 



20) 



Based on a technique from algebraic geometry called prime ideal decomposition and 
its implementation in the software SINGULAR [29], one can decompose the system (2.20) 
into its prime ideals which can be solved separately and give rise to disconnected branches 
of the general solution. Only one of these branches happens to contain dS extrema and 
these organise themselves into a one-parameter family of solutions. This dS line is plotted in 
figure 3 together with the value of the rj parameter as a function of the scan parameter, which 
shows the presence of tachyons at every point. On the other hand, this over-constrained 
situation already suggests that one might not be able to tune the masses of the six isotropic 
real scalars and get stable dS points. Further in figure 4 we show the normalised energy 7 
and its corresponding critical value saturating the sGoldstino bound in (2.14). This analysis 
manifests that the averaged sGoldstino (over real and imaginary parts) is never tachyonic 
all along the dS, with 7 that stays orders of magnitude away from the critical value and only 
approaches it in the singular regime in which 5 goes to infinity^. 



^For further discussion on the presence of tachyons in these models, and the 5 hmit in particular, see 
ref. [16]. 
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Figure 3: Left: the one-parameter branch of dS solutions of geometric type II A isotropic flux 
compactifications . The value of the energy Vq is given as a function of the scan parameter 
S = exp(Ai). The curve crosses the Minkowski line at 6 ^ 5.127. 

Right: the value of the rj parameter within the isotropic sector as a function of the scan 
parameter 6 already reveals the presence of tachyons all along the dS branch. 



-6 - 

-7 - 



Figure 4: Comparison between the value of'y along the dS line and 'y critical = 
maximal value of -y for which the averaged sGoldstino mass is not negative. 



3 



1, i.e. the 



3 Towards stable dS solutions 

In the previous section we have reviewed the case of geometric type IIA compactifications 
as an over-constrained example to which one can think of applying the method presented in 
section 2.2. Now we would like to move to the most interesting situation in which M exactly 
equals 2N . In our isotropic setup (A^ = 6) this implies the necessity of having 12 fluxes. 
By taking a look at tables 1 and 2, one realises immediately that this requires including 
non-geometric fluxes both in type IIA and in type IIB. 
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3.1 A type IIB setup with generalised fluxes 

The setup we focus on is the IIB duality frame with 03 and 07-planes with generalised 
fluxes defining a locally geometric background, i.e. F3, H, Q and P fluxes. The induced 
superpotential is the one given in (2.11). The aim of this section will be solving the equations 
of motion in the origin (2.18) for the isotropic scalars in order to find the most general set 
of extrema in this case and study the stability features thereof. 

We will analyse here the special sub-case of a generalised type IIB background with only 
F3, H and Q fluxes and vanishing P flux. The general superpotential in this case is given in 
(A.l). This choice is motivated by the fact that, in this non-geometric setup, a dS critical 
point was found in ref. [17] with the interesting property of being stable with respect to the 
six isotropic scalars. This supergravity solution, though, reveals the presence of tachyons 
in the non-isotropic directions of the scalar manifold. We will now solve the equations of 
motion for the scalars in this setup and show that the solutions span a parameter space of 
dimension six, in which the stable dS point of [17] sits. 

In general, given such a 6-dimensional parameter space, one can only hope to have control 
on the eigenvalues of the mass matrix in the six isotropic directions and there is no freedom 
left to stabilise the extra 8 non-isotropic scalars which could then remain tachyonic. 

The method for flnding the general exact solution to the field equations in this non- 
geometric type IIB setup is presented in appendix A. 

3.2 Statistical scan and tests of RMT 

In the previous section we have motivated our interest in the study of non-geometric type IIB 
compactifications with F3, H and Q fluxes. This seems to be the minimal setup for producing 
stable dS solutions at least for what concerns the mass spectrum within the isotropic sector. 
At least qualitatively we have argued that this situation in which N — 6 fields are retained 
can be regarded as a toy example for extrapolating some universal features of the probability 
distribution of stable critical points for potentials with N fields in general. 

Analysing the probability distribution of stationary points of different types {i.e. min- 
ima, maxima or saddle points) for classes of scalar potentials is equivalent to studying the 
distribution of the eigenvalues of its Hessian matrix evaluated at those given points. In the 
case in which one assumes the scalar potential to be a random function, its Hessian matrix 
will be a symmetric matrix with random entries. 

The distribution of the eigenvalues of a random matrix is the subject studied by Random 
Matrix Theory (RMT). This framework has been used in refs [30-33] in order to predict the 
heavy suppression of stable dS vacua in string-inspired supergravity models with N large 

13 



enough. 

Let us now go back to our general solutions presented in section 3 and appendix A. Even 
though this set of exact solutions is obtained analytically, one immediately realises that 
scanning a parameter space of dimension six in order to find stable vacua requires a different 
approach, since the analytical beahviour of the eigenvalues of the mass matrix as a function 
of these parameters is out of reach. 

What we did was generating a statistically significant sample of random values for the six 
supersymmetry breaking parameters {Aa, Ba}ci = sTU^ thus yielding an equally big sample 
of random points in the space of solutions. In figure 5 we show the resulting probability 
distribution of tachyon-free critical points as a function of the scanning parameter^ 7 rep- 
resenting the ratio between the positive contribution to the cosmological constant and the 
AdS scale. 

This result shows the presence of a peak of tachyon-free points in the AdS region with 
small 7 but away from the supersymmetric point at 7 = 0. This could be explained by 
observing that the mass spectrum in the supersymmetric case is completely determined by 
the masses m^^ of the chiral fermionic superpartners of the scalars. When these are massless, 
the corresponding partener scalars sit at the conformal value — Non-zero fermionic 
masses will produce linear corrections to the off-diagonal elements of the scalar mass matrix 
and quadratic ones to its diagonal entries. By studying in detail the sign of the eigenvalues, 
one finds the presence of at least one tachyon (though obviously still above the BF bound) 
as a generic feature, whereas introducing supersymmetry breaking makes it slightly easier 
to arrange for all positive eigenvalues. 

Furthermore, there appears a second sharp peak in figure 5 in the AdS region but with 
values of 7 close to 1. We do not have a clear explanation for this at the moment, though it 
seems to indicate that there might be a natural way of saturating the sGoldstino bound for 
the 77 parameter in the AdS region for a certain value of 7. 

Finally, we observe a tiny tail of stable dS points with the very specific feature of having 
7 slightly above 1. This seems to be in hne with the results in refs [34,35] where the relation 
between flux-induced potentials and RMT predictions is discussed, the final outcome being 
that string compactifications strongly suppress stable dS vacua with large values of the 
cosmological constant. 

In what follows, we would also like to make a comparison between our present analysis 
and the RMT-based results. It is important to keep in mind that our supergravity models 
describing a background with generalised fluxes only allows for a small number of random 

^Please not that 7 = 7 + 1,7 being the normaUsed energy parameter introduced in the context of the 
sGoldstino analysis. 
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parameters that equals the number of diagonal entries of the mass matrix, whereas all the 
remaining off-diagonal ones are determined in terms of these by the equations of motion. In 
other words, while the RMT is dealing with a more or less arbitrary random superpotential, 
we are considering a constrained string inspired subset. 

We do not expect the RMT description to be applicable when the supersymmetry break- 
ing scale and the AdS scale are of comparable size (7 > 1). However, it is interesting to 
make a comparison in the far-right region where 7 is significantly larger than the Minkowski 
value of 1. Using the estimate [32,36] 

r log3 , log(2V3-3) logA^ 

V = exp — + '- N 0.0172 

^ 4 2 24 

one finds an estimated fraction 4.36 x 10""^ of stable dS for N — 6. This would predict 
roughly 25 to 30 stable solutions out of the total of ~ 6.3 x 10^ dS solutions with 7 > 1.5 
in our scan. In contrast, we were unable to find a single example of stable dS in this region. 
This significant deviation from the RMT prediction confirms our expectation that there is 
actually no regime in which our Af — 1 supergravity potentials, inspired by non-geometric 
string backgrounds, fit an RMT approximation. 

If one wants to reproduce the RMT predictions for stable dS vacua in the high 7 region, 
it is clear that one needs to include many more random parameters. The stringy origin of 
these remains, however, an open issue. 

As a motivation towards further work in this direction, we would like to give a first look 
at the N — 14: case corresponding to studying the full SL(2)^ model relaxing isotropy. The 
analgous setup allowing for enough freedom in parameter space will now have 28 fluxes. 
By taking a look at the last column of table 1, one observes that also in this non- isotropic 
situation one needs to consider non-geometric backgrounds in any duality frame. As an 
example of this, we chose a type IIB background with the most general gauge fluxes (F3 
and H) plus those Q flux components appearing in superpotential terms with up to linear 
dependence on C/j (The first two types of Q flux in table 1). These fluxes turn out to be 
exactly 28 in total. 

In figure 6 we give an idea of the results for the random scan in this case. These pre- 
liminary results were obtained by means of a limited sample of data and we leave a detailed 
analysis for future purposes. Such a study will certainly be needed in order to investigate 
the presence of fully stable dS vacua, if any. However, the present diagram of the proba- 
bility distribution of tachyon-free points w.r.t. all real moduli roughly shows some common 
features with the N — 6 case in the AdS region. 



(3.1) 
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Figure 5: Distribution of the fraction of tachyon-free critical points of the isotropic scalar 
potential with N = 6 fields as a function of the uplift parameter 7 = ^|p^p- ■ Data produced 
from the analysis of 10'^ random points. The zoomed diagram in the upper-right corner shows 
the presence of a stable dS tail with 1 < 7 ^ 1.02. This tail was analysed using a sample 
bigger by a factor of 10. 



3.3 Analysis of the stable dS region 

Now we would like to carry out a complementary analysis w.r.t. the statistical scan presented 
in the previous section. The aim is that of charting the region of stable dS inside parameter 
space in the JIB case previously introduced and studied. In this way, one can think of 
explicitely construct the concrete form of the conceptual picture given in figure 2 which 
applies in this case. 

Since the parameter space of solutions has dimension six, we studied different two- 
dimensional slices thereof. We have fixed four out of six supersymmetry breaking parameters 
to the values corresponding to the stable dS solution in ref. [17] and moved about in the 
remaining two-dimensional slice. We have constructed the level curves of the scalar po- 
tential around such a point and the level curves of the minimal eigenvalue of the mass 
matrix. Finally, we have studied the overlap between the dS region and the tachyon-free 
region thus locating stable dS points in parameter space. The results of such an analysis 
are collected in figures 7 and 8, which represent the {^{DsW),Q{DsW))-pla.ne and the 
{^{DsW), ^{DTW))-p\ane respectively. 
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Figure 6: Distribution of the fraction of tachyon-free critical points of the full scalar potential 
with N = 14 fields as a function of the uplift parameter 7 = ^|^^p-- Data produced from 
the analysis of 5 x 10^ random points. Although the statistical significance is not sufficient 
for detecting any stable Minkowski or dS tail, the general features in the AdS region are 
compatible with the results obtained for N = 6. 
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Figure 7: The parameter space of solutions projected on the {Ai, Bi) plane. Left: Level curves 
of the cosmological constant; the lower-right corner filled with lighter colours and delimited 
by the Minkowski (blue dashed) line corresponds to the dS region. Middle: Level curves of 
the 7] parameter; the tachyon-free region corresponds to the upper-left part filled with lighter 
colours. Right: The tiny region of overlap in parameter space corresponding with stable dS 
is zoomed in here. 

The striking pecuharity of the dS region and the tachyon-free region is that they happen 
to be almost each other's negative and this also holds for all the other two-dimensional slices 
that are not given here. This always makes the corresponding overlap region look like a very 
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Figure 8: The parameter space of solutions projected on the {Ai, A2) plane. Left: Level curves 
of the cosmological constant; the dS regions are the ones filled with lighter colours below by 
the Minkowski (blue dashed) lines. Middle: Level curves of the rj parameter; the tachyon-free 
regions are filled with lighter colours. Right: The tiny region of overlap in parameter space 
corresponding with stable dS is zoomed in here. 

narrow band. By moving away from the stable dS point in which the plots are centered, the 
stable dS band tends to persist. This may be regarded as an indication for the stable dS 
region to have the shape of a thick sheet in six dimensions, in which one of the dimensions 
is much narrower than the other ones. Studying the topology of such a stable dS region 
and in particular, as a first thing, understanding whether it consists of or more connected 
components, requires some further detailed investigation that we leave for future work. 

Furthermore, this complementary analysis corroborates what the random analysis was 
already suggesting, i.e. that the very tiny stable dS regions always occurs very close a 
Minkowski critical point and hence with values of 7 slightly above 1. 

4 Conclusions 

In this paper we have offered a framework for systematically analysing the space of solutions 
of type II string compactifications compatible with A/" = 1 supersymmetry. The supergravity 
effective description realises string dualities as global symmetries, and therefore it comprises 
different descriptions in different duality frames in a single universal formulation. 

Supersymmetry breaking parameters allow one to predict the dimensionality of the space 
of solutions by means of a very simple counting argument. In this context we were able to 
revisit the case of geometric type IIA compactifications. We find that the space of solutions 
does not seem to contain enough degrees of freedom for getting rid of the tachyons even 
within the isotropic sector. 
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Subsequently, we applied this novel prescription to investigate less constrained situations 
such as isotropic IIB compactifications with generalised fluxes. A statistical scan as well 
as a charting of two-dimensional slices of the parameter space reveal the presence of thin 
sheets of stable dS regions with small cosmological constant (compared to the AdS scale). 
Generically, since our parameter space is just big enough for controlling the eigenvalues of 
the mass matrix in the isotropic sector, we expect that the remaining non-isotropic sector 
will contain tachyons. This we have checked in some explicit cases, even though we have not 
performed a full analysis. 

However, by considering non-isotropic setups with a larger parameter space, our general 
expectation is that there will exist a non-empty stable dS region. This region will probably be 
considerably thinner than its isotropic counterpart. Indeed, a first random scan of the non- 
isotropic case (although not with enough statistical significance), has exhibited no examples 
of tachyon-free dS critical points. This means that the issue of the existence of a fully stable 
dS vacuum, w.r.t. to all fourteen non-isotropic moduli, still remains to be addressed. We 
hope to come back to this point in the future. 

Our work shows that stable de Sitter vacua are much less frequent than one naively 
would expect. In particular, we have observed the necessity to add at least some flux that 
is non-geometric from both the type IIA and the type IIB point of view. Furthermore, 
our analysis hints at tantalizing conspiracies, which suggest that dS vacua, if at all present, 
have unnaturally small values (compared to the AdS scale). It would be interesting to 
analytically understand the reason for these conspiracies, and to investigate their possible 
phenomenological consequences. 
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A All the solutions of type IIB with F, H and Q 



The general form of the superpotential is obtained by taking (2.8) and switching of all the 
couplings but those ones in Pp, Ph and Pq. This yields 



(A.l) 



+ 3r(co + Ci[/-C2t/'-C3t/3) 

After solving the "broken" supersymmetry conditions (2.15) in terms of {62, &3, Cq, Ci, C2, C3} 



via 



&2 
&3 

Co 
Cl 
C2 

C3 



I (3ai - aa + 26o + ^1 - A2) 
I (-ao + 3a2 + 6&1 - Si + S2) 
(9ai + 03 + 46o - ^1 - ^2 - 2^3) 
(9ao - 1502 - 12&1 - 3Si - 3^2 - 2S3) 
^ (9ai - 7a3 - 46o - 5^1 - ^2 + 2A3) 
3ao - 3a2 + 126i - 2>Bi + B2 + 2B3) 
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(A.2) 



one finds 



(65i + 3^2 - 3^3)00 + 3(6Ai + 3^2 - ^3)^1 - 9(2Bi + B2 - ^3)02+ 
-(6A1 + 3^2 - ^3)03 - 8A360 + 24^3^1 - QAj - Ai{3A2 + As) - 6Bl+ 
-W1B2 - WiB-i -Al + A2A3 - Bl + 352^3 + 2Ai + 2Bl = , 

3(2Ai -A2 + A3)ao - 3(6fii - 3^2 + ^3)01 - 9(2Ai - A2 + ^3)02 - Sfig&o - 24A361+ 
+(6Si - 3^2 + ^3)03 - Ai(3S2 + Bs) + 3Si^2 + 35^3 + ^2^3 - 352^3 = , 

3(3Si + 4^2 + 3^3)00 + 9(3^1 + 4^2 + ^3)01 - 3(9Si + 12^2 + ^3)03+ 
-(9^1 + 12^2 - 5^3)03 + 8A3&0 - 24S361 + %Al + Ai(-15A2 + ^3) + 9^2+ 
-I5S1S2 + 3S1S3 - AAl - A2AS - A.BI - 3S2S3 + 2Al + 2S| = , 

-9(Ai + ^3)00 + 9(3Si + 53)0,1 + 3(9Ai + As)a2 - (9fii - 5^3)03+ 
+8E360 + 24A361 + Ai(352 + Bs) - 3B1A2 ~ 3B1A3 - A2B3 + 352^3 = , 

-9(5i - ^2 - 2fi3)ao + 3(-3yli + 3^2 + 2^3)01 + 3(95i ~ B2 - 2^3)02+ 
+(3Ai + 5^2 - 18^3)03 + 8A260 - 24fi2&i - 3^2 + Ai(-246o + 4^2) - 9^2+ 
+72Bibi + I2B1B2 -Al- AA2AS - Wl - AB2B3 + AAl - 4S| = , 

3(3^1 - 3^2 - 2^3)00 - 3(3Si - 3^2 + 10^3)01 - 3(9^1 -A2 + 10^3)02+ 
+(3Si + 5^2 - 6^3)03 + 8S260 + 24A261 - ^i(-6Si + 726i + 8^2)+ 
-24Si6o + 2^252 + 4A2S3 - 4S2A3 + 8A3S3 = . 
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The above set of equations can be then solved with respect to {oq, Oi, 02, 03, bg, bi} which 
only appear linearly. This procedure combined with (A. 2) gives rise to to the general solution 
expressing the 12 fluxes as a function of the six supersymmetry breaking parameters. 
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